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We show how to engineer an extensive range of non-linear Hamiltonians for a nano-mechanical 
resonator. The technique requires only a time dependent drive applied to a Cooper-pair box or 
second oscillator to which the nano-resonator is coupled. This method allows one to generate a large 
number of non-classical states, as well as Hamiltonians whose classical counterparts are chaotic. 
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Nano-mechanical resonators can now be built with fre- 
quencies of 100MHz and quality factors of 10 5 [ll. Cou- 
pled with recent advances in the measurement [2j and 
cooling 0, 0, 0, 0, 0] of these systems, this has opened up 
the exciting prospect of observing quantum behavior in 
mesoscopic mechanical systems for the first time. To fully 
explore the quantum behavior of nano-rcsonators, and 
consider exploiting this behavior in future applications, 
one must prepare these systems in highly non-classical 
states, which in turn requires the action of a non-linear 
Hamiltonian. A number of approaches have been taken 
to this problem. One is to engineer non-linear couplings 
to @], or engineer a reservoir via 0], a Cooper-pair Box 
(CPB), but so far only a very limited set of Hamiltoni- 
ans have been obtained in this way. A second approach 
is to use measurement and/or feedback 11 [• While 
this has potential, low-noise measurements are, at least 
presently, much more difficult to perform than Hamil- 
tonian control. A third more exotic approach borrows 
from quantum optics, and involves embedding a quan- 
tum dot on the resonator, which in turn interacts with 
laser fields Here we show that a large range of non- 
linear Hamiltonians, and a correspondingly large range 
of non-classical states, can be engineered by coupling the 



resonator to a single CPB 0, m m or second nano- 
resonator. This does not require dissipative, measure- 
ment or feedback processes, but simply the application 
of time-dependent control fields to the auxiliary system, 
which is a relatively straightforward process. Note that 
our purpose here is to engineer nonclassical states of the 
resonator alone, and not an entangled state of the cou- 
pled systems. The latter can be achieved with a simple 

m 



linear Hamiltonian, as shown in 



1( 



This work was inspired by a 2004 article of Lloyd, Lan- 
dahl and Slotine [17] . in which they point out that, at 
least in principle, a single interaction of any system with 
a second simple system should suffice to create any dy- 
namics between the two systems, so long as the interac- 
tion operator, G, docs not commute with the system's 
Hamiltonian, H . This is because commutators of G and 
H will in general generate a complex algebra. These au- 
thors did not give a method for achieving any specific 
evolution, however, and there is probably no general an- 
swer. For the problem of state-preparation, the task is 



complicated by the fact that the interaction will in gen- 
eral entangle the two systems. To prepare a state of the 
primary system, we must entangle and then disentangle 
them in such a way as to leave a non-trivial operation on 
the primary system. One is reminded of attempting to 
flip an edge piece in Rubies cube while leaving the other 
facets intact. 

Here we show how to generate a significant range of 
non-linear Hamiltonians for a nano resonator by perform- 
ing time-dependent rotations on a single qubit to which 
it is coupled. We note that we are not the first to show 
that general quantum states of an oscillator can be en- 
gineered by time dependent driving of a low-dimensional 
auxiliary system. A scheme to do this was first proposed 
by Law and Eberly (l8j |. and another more recently by 
Santos These methods should also be realizable in 
nano-resonators, although they are very different from 
the technique presented here; these schemes are most ap- 
propriate for creating arbitrary states with small num- 
bers of phonons, since they become increasingly com- 
plex as the number of phonons increases. Conversely, 
the method described here is most appropriate for cre- 
ating states of arbitrarily large phonon number, such as 
mesoscopic Scrodinger-cat states, and states generated 
by chaotic systems such as the Duffing oscillator [io| . 
Further, our approach is not limited to two-level auxil- 
iary systems. We show how to generate many non-linear 
Hamiltonians with an auxiliary oscillator, which may well 
be important in nano-mechanical applications, and the 
method could certainly be extended to other systems. 

We consider first a Cooper-pair Box as the auxiliary 
system. If we electrically bias the nano-resonator and 
place it adjacent to a CPB, then the two systems in- 
teract via the Hamiltonian H\ = axa z , where a is the 
coupling constant determined by the geometry and bias 
voltage, x = a + is the position of the resonator and 
a z is the Pauli operator for the CPB qubit in the charge 
basis. If we wish to move into the interaction picture 
to eliminate the free rotation of the resonator, we can 
preserve this time independent interaction by modulat- 
ing the coupling strength at the resonator frequency and 
making the rotating wave approximation. If instead we 
far-detune the CPB from the resonator, and operate the 
CPB at the charge degeneracy point, then the interaction 
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becomes H 2 — a'Na x where N = a^a [21|. We can ob- 
tain an effective oscillator Hamiltonian proportional to N 
either by changing the modulation frequency of the cou- 
pling strength, or by employing H 2 with the qubit in a a x 
eigenstate. Finally, we can perform arbitrary rotations of 
the CPB qubit by changing an adjacent gate voltage 
These are the basic interactions and operations that we 
wish to use to create non-linear Hamiltonians for the res- 
onator. We will do this by using the interactions, along 
with sequences of rapid rotations applied to the qubit. 
In what follows we set h = 1 to simplify the notation. 

The first tool we employ is the following: given the in- 
teraction Hamiltonian Acjj , where A is an operator of the 
resonator and aj is a Pauli operator for the qubit, one can 
generate the evolution e~ lAnrT . To do so one first applies 
a rapid rotation to the qubit given by e (where m is 
defined by [m X n] • a = Uj ) , allows the systems to inter- 
act, and then applies the inverse rotation. The result is 
the desired unitary, since g-'^'^e^ = e -iAn <r 

Our second tool is the Zasscnhaus formula, which may 
be regarded as complementary to the more well-known 



BCH-formula 22, 23 |, and allows one to write an expo- 
nential of the sum of two operators as an infinite product 
of exponentials of the operators and their commutators: 
P (A+B) _ p A p B p c 2p c 3 .. . ^| xhe second and third 



order terms in the product expansion are 



C 2 = -(1/2) [A, B], 
C 3 = (1/3)[B, [A, B}] 



(1/6) [A, [A, B] 



(1) 
(2) 



and further terms are given by a recursion relation [25[ . 

We now assume (without loss of generality) that the 
interaction between the resonator and auxiliary qubit is 
aBa x , and use the above procedure, in tandem with 
the Zassenhaus formula, to generate an evolution pro- 
portional to the square of B. To do so we first change 
the qubit representation to <r = [a , x ,a y ,d z ) = ([a x + 
a z }/\/2, a y , [<r x — a z ]/\/2). The evolution due to the in- 
teraction during a time interval At becomes exp(iA(a x + 
<7 2 )), were we have defined A = aAtB to simplify the no- 
tation. We then expand this expression to second order 
using the Zasscnhaus formula, and multiply on the left 
by e -iA&* e -iA<? x to 0D t a in 



AAo 



_ e -iA A a y /2 



(3) 



valid to second order in A. This gives us the following 
procedure for generating the Hamiltonian A 2 . We first 
apply control rotations to the qubit to prepare it in an 
eigenstate of a y . We then wait a sufficiently short time 
interval during which the systems interact via aA<j x . To 
complete the operation we apply the unitaries e _4j4crx and 
e -iA(Jz ^ ky usm g the procedure described above. During 
this sequence the two systems are entangled and then 
disentangled, and the resonator undergoes the non-trivial 
evolution e ± l ( aAt ) B . Here the sign is given by the choice 
of the initial qubit eigenstate, and the evolution time 



must be small compared to the eigenvalues of aB. Ap- 
plying this full sequence repeatedly then generates this 
evolution for an arbitrary time t. 

Note that using the above procedure to generate the 
evolution A 2 requires that the time interval in each repe- 
tition be sufficiently small that the third and higher order 
terms in the Zassenhaus expansion can be ignored. We 
can therefore increase dramatically the procedure's effi- 
ciency by obtaining an expression similar to Eq.([3|) that 
is valid to third order instead of second. We can do this 
by noting that if we flip the signs of both A and B in the 
Zassenhaus formula, the sign of the second order term 
does not change, but that of the third order term does. 
We can therefore cancel the third order term by applying 
the sequence again with the signs flipped. This gives 



e iA ~°*e iAa ~*e 



iA °*e 



%Aa x —iA(f x iA<r x 



-iA 2 a y /2 



(4) 



which is now valid to third order in A. 

In creating the effective Hamiltonian A 2 , we have also 
obtained the effective interaction A 2 a y . We can there- 
fore apply the above procedure to this new interaction 
to obtain the effective Hamiltonian A 4 . In theory this 
procedure allows us to generate all even powers of A. 

To create odd powers of an operator A, we can use 
the same method, but this time expand the Zasscnhaus 
formula to third order. This gives 



^iA <7 y jl piAtTz giAcrx ^—iAa x 



-iA 3 [25- x -(fz]/12 



(5) 



To apply the first unitary on the left hand side, we can 
use essentially Eq.Q, but we must change the sequence 
on the left hand side a little so as to flip the sign on the 
right hand side. The operation we require is 



giASz — iAdx e iAa z 



AA 2 d y /2 -iA 3 [2a x +a z ]/12 



Combining Eq.([5]) and Eq.© we obtain finally 

piAoz g— iAa x piAcr z piAa z piAtfx g— iAa x _ g — iA 3 <f z /6 



(6) 



(7) 



which is an explicit pulsc-scqunce for generating A 3 . We 
can also increase the efficiency of this procedure, as we 
did for the A 2 sequence above, by expanding the Zassen- 
haus formula to forth order, and choosing a sequence that 
eliminates the forth order terms. We now have explicit 
methods to generate second, third and forth powers of 
any interaction operator, and it is clear that by expand- 
ing the Zassenhaus formula to higher orders, we can gen- 
erate higher powers if desired. 

As discussed above, we have at our disposal interac- 
tion Hamiltonians with A oc x and A oc N. By using 
both together we can create all real linear combinations 
of powers of x and N. Note that the squeezing Hamil- 
tonian is in this set, being a 2 + a^ 2 = x 2 — 2N. To ap- 
ply both interactions together (or sequentially, as these 
are equivalent to first order) to obtain the desired lin- 
ear combinations, the simplest method would probably 
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FIG. 1: Wigner functions for various states engineered with the techniques discussed. Plot key: luminosity is proportional to 
absolute value, with blue positive and green negative, (a) A squeezed Schrodinger-cat state; (b) three superposed coherent 
states; (c) two superposed displaced number states; (d) a state generated by the Duffing oscillator after ten drive periods. 



be to have the resonator interact with two CPB's, one 
on-resonance and the other detuned. 

So far we have only used time-dependent control per- 
formed on the qubit. It is worth noting that if we also 
apply rapid "rotations" to the resonator, given by e t6N , 
then we can generate an effective interaction proportional 
to a linear combination of x and p by using the relation 

e -ieN e i\x<r, e i8N = e i[cos(e)x+sin(e)p}a z _ This aUows us to 

further extend the accessible set of non-linear Hamilto- 
nians. Rapid rotations of the resonator are equivalent to 
shifting the resonator frequency, which could be achieved 
by using the interaction XNa z , and modulating A. 

We now consider using an oscillator as the auxiliary 
system. Coupling to an auxiliary nano-resonator via 
XxX is straightforward [2(|, where X = (b + tf) is the 
position of this resonator. Using a superconducting L-C 
oscillator as the auxiliary system one can obtain the inter- 
action fiNb^b We allow ourselves two simple control 
operations on the auxiliary oscillator. The first is a linear 
drive, and the second is a frequency shift (phase-space ro- 
tation). Given an interaction aA(b + b^), where as before 
A is a resonator operator, the rotations allow us to obtain 
all effective interactions of the form a.A{be~ l + b^e -10 ). 
The Zasscnhaus formula provides the exact relation 



iAX iAP -iA(X+P) 

e e e y ' 



(8) 



giving a pulse sequence to generate A 2 . Because this 
relation is exact, the pulse sequence need only be applied 
once to obtain A 2 , in contrast to the sequences above. 
This advantage follows because the Zassenhaus sequence 
C n terminates for X and P. As a result, we cannot use 
X and P to obtain higher powers of A. 

We can obtain all even powers of N by using 
e -i\Ntfb+iab+ia*bi and expanding as before with the 
Zassenhaus formula. Odd powers can be obtained by 



sandwiching e 



iNb^'b 



between linear driving pulses e ±zaX 



to give e lNb b+laNP ; and once again applying Zassen- 
haus. As an example, the sequence we obtain to gener- 
ate N 3 (valid to third order) is U{-X)U(X) = e l ( 2 /3)A 2 w 3 
where 



U(X) 



JXNP iNb*b -iN(b*b+\P) 



We can use the above methods to engineer a range of 
non-classical states, including: 

Squeezed Schrodinger-cat states: A "Schrodinger Cat" 
is a superposition of two different coherent states. To 
generate a superposition of two squeezed states, we first 
cool the resonator to the ground state, and then apply 
the squeezing Hamiltonian a 2 + a <f2 to produce a squeezed 
vacuum state. Next we apply a classical linear drive, 
given by the Hamiltonian eta + a*a^ to shift the state in 
phase space. Finally we apply the Hamiltonian A(aJa) 2 
for a time t = n/(2X), and this generates a superposition 
of two squeezed states symmetrically placed about the 
origin [28|, [29[ . The Wigner function for such a state is 
displayed in Fig. 1 (a). 

Multi-cat states: To generate a superposition of more 
than two coherent states we can again use the Hamilto- 
nian A(a^a) 2 . Applying this Hamiltonian to an initial co- 
herent state for times t = 7r/(3A) and t = 7r/(4A) creates 
superpositions of three and four coherent states respec- 
tively. The Wigner function for the former is shown in 
Fig. 1 (b). 

Displaced-number cat states: A displaced number 
state is a number state that has been shifted in phase 
space [Io[ • We can create number states in a straightfor- 
ward way by starting with the ground state, and using 
the qubit to successively add single phonons to the res- 
onator. To do this we employ the interaction Xxo~ x , which 
under the rotating wave approximation simply generates 
an exchange of excitation between the qubit and the res- 
onator. If we start the qubit in its "upper state" |1), and 
the resonator in the ground state, |0), then after a time 
7r/A the resonator is in state Repeating this one can 
obtain any \n) . We can then displace these number states 
by applying a coherent drive, and subsequently form a 
superposition of two or more by applying A(ata) 2 . A su- 
perposition of two displaced number states with n = 2 is 
shown in Fig. 1 (c). 

Preparing non-classical states is not the only moti- 
vation for generating non-linear Hamiltonians in nano- 
electro-mechanical systems. Such Hamiltonians also al- 
low us to engineer specific kinds of dynamics for these 
systems, including those whose classical counterparts are 
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chaotic. These dynamical systems are of special interest 
from the point of view of the quantum-to-classical tran- 
sition, and especially the emergence of classical chaos. 
This is of particular interest for nano-resonators, because 
1) it has been shown that the continuous measurement of 
the position of a non-linear resonator is sufficient to in- 
duce the transition from quantum dynamics to classical 
chaos HEll, and 2) a continuous position measurement 
of a nano-resonator has already been realized experimen- 
tally 0. An example of a chaotic system with a single 
degree of freedom is provided by the Duffing oscillator, 
whose Hamiltonian is aa)a — Xx 2 + /ix 4 + A cos(wi)x [2(| . 
This Hamiltonian can be engineered using the above tech- 
niques, since we can generate a' a by detuning, x 2 and x A 
using the methods above, and the final term is merely a 
time dependent linear drive. As an example, we illustrate 
a typical state generated by the Duffing oscillator after 
10 periods of the drive, from an initial coherent state, 
and using the parameters in [2oT | with h = 1. 

We now discuss realizing the above techniques with 
current experimental parameters. Our goal is to gener- 
ate dynamics that is significantly faster than the damp- 
ing rate of a nano-mcchanical resonator, which is cur- 
rently about T = 10 4 s~ 4 . The control operations are 
applied to the auxiliary system (CPB or auxiliary oscil- 
lator) by switching the voltage on an adjacent gate. It 
is not difficult to achieve switching times on the order 
of a nano-second, and significantly faster switching can 
be achieved with advanced electronics. Since the dynam- 
ics engineered in the resonator must be smaller than the 
switching rate, this is a limiting timescale. 

Consider as a concrete example preparing the cat state 
(|a) + \-a))/y/2, with (a\N\a) = 10. To do this using 
an L-C oscillator we require the interaction /iNb'b, and 
in this case fi can realistically be as high as 10 8 s _1 [27j |. 
(For a CPB, the equivalent interaction strength can be at 
least 10 6 s~ 1 [n|.) If we set fj, = 2.5 X 10V\ and gener- 
ate the Hamiltonian using 50 third-order pulse sequences, 
so as to limit unwanted higher-order terms to 2%, we can 
prepare the state in a time of r = \/50 x 27r/(4^) = 5/is. 
A switching time of Ins should be adequate for this pur- 
pose, as the whole procedure will require a few hundred 
pulses. (With faster electronics we could use more pulses 
and achieve higher accuracy.) To gauge the effect of de- 
coherence, we now simulate numerically the action of the 
N 2 Hamiltonian with the addition of thermal damping at 
r = 10 4 s _1 and a temperature of lOmK. For a 100MHz 
resonator we find that the damping reduces the height 
of the interference fringes from that of a perfect cat by 
about 50%. This would therefore be an appropriate sce- 
nario for probing environmental decoherence. The L-C 
oscillator is also subject to decoherence at a similar rate, 
but we expect this to have less effect on the final state, 
since it is repeatedly disentangled from the resonator dur- 
ing the pulse sequence. On the other hand, decoherence 
times for CPB's arc considerably shorter. While these 



have been steadily improving, they currently stand at 
about 1 fis [32j], and this might well interfere with the 
preparation process. A detailed investigation of the ef- 
fects of decoherence in the auxiliary system will be the 
subject of future work. 
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